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Abstract 

The construction of a generalized (higher-order) nonhnear thermo- hydrodynamics, based on a 
nonequihbrium ensemble formalism has been presented in the preceding article. The working of 
such theory is illustrated in the present one. We consider here the case of two ideal classical fluids 
in interaction between them, and the nonlinear equations of evolution for the density (a hyperbolic 
one) and for the velocity field describing motion under fiow are derived. Also, now at the quantum 
level, it is described the nonlinear transport in the fluid of electrons in doped polar semiconductors, 
and comparison with Monte Carlo calculations and with experimental data is done, obtaining very 
good agreement. 
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I. INTRODUCTION 

In the preceding article (hereafter called as I, and when equations present in it are in- 
dicated in what follows, we write the number preceded by I) 1], it was described the con- 
struction of a generalized (higher-order) and nonlinear thermo-hydrodynamics, based on 
mechano-statistical foundations. In that way it was established a kind of unification of the 
microscopic dynamics and a mesoscopic hydrodynamics. 

In this follow-up article we attempt to illustrate the functioning of the theory by applying 
it to a couple of somewhat simplified situations. Both cases are such that we can describe 
their underlying mechanics in terms of an individuals single-particle approximation. They 
are: 1. The case of two ideal classical fluids in interaction between them, one being the 
main system of interest driven away from equilibrium, and the other acting as a thermal 
reservoir; the associated nonlinear generalized hydrodynamics, however restricted to a study 
of a first-order one, is derived in section II. 2. A similar case, but now composed of two 
quantum gases, namely, conduction band electrons in a doped polar semiconductor under 
the action of an eletric field, and the phonon system acting as a thermal bath, is dealt with 
in section IV. Moreover, the study of section II is complemented in section III with the 
analysis of diffusion-advection motion governed by a nonlinear hydrodynamic equation. 

II. TWO CLASSICAL FLUIDS IN INTERACTION 

Let us consider two ideal classical fluids with interaction between them, with the Hamil- 
tonian given by 



m being the mass of the particles in the system of interest, M that in the second system 
taking as a thermal bath at temperature To, and the interaction being of the type of central 
forces. Vext stands for the interaction energy with an external source, and it is to be under- 
stood that the bath is constantly in equilibrium with an external ideal reservoir. We analize 
the hydrodynamics of this system in a first order approximation, meaning to take as basic 
variables of the system of interest the kinetic energy, the density, and the first flux (current) 




(1) 
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of the latter, with the associated dynamical quantities being 



Mr) = E|^^(r-r.), (2) 



n 



r) = J2S{r-r,), (3) 



I'^W = E^^(r-r^)' (4) 

fit 

j 

with the understanding that all three are defined in F-phase space. The thermal bath is 
assumed to remain constantly in equilibrium with a reservoir at temperature Tq, and the 
system subjected to external forces driving it out of equilibrium. 

According to I (the previous article) the auxiliary ("instantaneously frozen") nonequilib- 
rium statistical operator is given by [cf. Eq.(1.6)] 

n{t, 0) - p{t, o)xpR, (5) 

where pji is the canonical distribution at temperature Tq of the thermal bath, and 

p{t, 0) = exp I -(f){t) - J dh [F^(r, t) h{v) + A{r, t) n{r) + F„(r, t) • in(r)] | , (6) 

where, we recall, (j){t) ensures the normalization and we have introduced the set of nonequi- 
librium thermodynamic variables (Lagrange multipliers in the variational derivation of the 
formalism), namely 

{Fh{r,t),A{r,t),F„{r,t)} , (7) 

and we also recall that the noncquilibrium statistical operator is given, in terms of the auxil- 
iary one of Eq.(5), by Eq.(1.4). The corresponding set of basic hydrodynamic macrovariables 
consists of 

h{r,t) = J dTh{r)p{t,Q), (8) 
n(M)^/dr.Mp(,0), (9) 



In{r,t) = J drin{r,t)p{t,0), 
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(10) 



where the integration is over the phase space of the system [cf. Eq.(1.31)]. These equations, 
Eqs.(8) to (10), are the nonequilibrium thermodynamic equations of state, meaning that 
they relate the basic macrovariables on the left, with the nonequilibrium thermodynamic 
variables of Eq.(7) present on the right-hand side in the distribution p of Eq.(6). 
Equations (8) to (10) take the form 

h{r,t) = J d'p^f{T,p;t), (11) 

n(r,t)= J rf^p/(r,p;t), (12) 

I„(r,i)= / d'p^f(v,p;t), (13) 
J rn 



where 



/(r, P; = ^ exp I A(r, t) - F,(r, t)^- F„(r, t) ■ ^ , (14) 



is a Boltzmann-style one-particle distribution function, and 



Z{t) = J d\ J d3pexp|A(r,t)-F;,(r,t)^ + F„(r,t)-^| . (15) 
Performing the calculations we find that 

where we have introduced the barycentric velocity v(r, i) [cf. Eq.(20) below], through the 
definition (see also Eq.(20) below) 

Fn(r,t)^Fh(r,t)mv(r,t). (17) 

Moreover, we do find that 

h{r,t) = u{r,t) + n{r,t)^mv^{r,t) , (18) 
for the energy, where u is the internal energy 
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introducing, via the definition = fc^ T*, the field of quasitemperature j2,y,yi or nonequi- 
librium kinetic temperature j^] T*{r,t), with the last term on the right of Eq.(18) being the 
drift-kinetic energy, and 

I„(r,t) = n(r,t)v(r,t). (20) 

It can be noticed that if n is the concentration of particles, it needs be satisfied the 
constraint that 

n^jSruir,t), (21) 

with n{r,t) given by Eq.(12). 

Let now go over the generalized hydrodynamic equations, namely those of Eqs.(1.35), 
(1.36) and (1.38) of I (previous article). We do have that 

|^Mr,t) + V-I„(r,t) = A(r,t), (22) 

where Ih is the fiux of energy and Jh is the collision integral; 

^n(r,t) + V-I„(r,t) = 0, (23) 

which is the conservation equation for the density; 

^I„(r, t) + V ■ If (r, t) = J„(r, t) + Fe.t(r, t) , (24) 

where In (r, t) is the fiux of fiux (second-order fiux) of particles, namely 

/^](r,t)= / rfr5^[P^P^j5(r-r,)p(t,0), (25) 

j 

with [pj pj] standing for tensorial product of vectors, rendering a second-rank tensor, 
Fe^t(r, t) stands for external forces [arising out of Vext in Eq-(l)] applied on the system driven 
it out of equilibrium, and Jn is the collision integral accounting for relaxation processes (to- 
wards the thermal bath): See Appendix A. The latter is composed of three contributions: 
two of them, one introducing a term proportional to the gradient of concentration (and 
responsible for a perturbational modification of the diffusion effect; cf. Ref.[6]) and another 
introducing nonlocal effects (correlations in space), are neglected and we conserve the main 
contribution given by 

J„(r,t) = -[e[2l(r,t)]-'l„(r,t), (26) 
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introducing the inverse of a second-rank tensor, 0„ having dimension of time and playing 
the role of a tensorial-in-character Maxwell-relaxation time 0,Q)S|, its explicit expression 
given in Eq.(A3) of Appendix A. 

Finally, to close the system of Eqs.(22) to (24) coupled to the set of Eqs.(16), (18) and 
(20) - the equations of state - we need to express the second-order flux of Eq.(25) in terms 
of the basic variables, and we obtain that 

mJPl(r,t) = n(r,t) [A;^ r*(r, t) I'^l + mn(r, t) [v(r, t)v(r, t)]] . (27) 

This second-order flux is related to the field of pressure tensor by 

Pl2](r,t) = m/|l(r,t) -n(r,t) [v(r, t)v(r, t)] = n{r,t) keT* {r,t) fl^^ = -u{r,t) l'^] , (28) 

3 

after using Eq.(19), recovering the local and instantaneous usual form of the hydrostatic 
pressure. 

This first-order hydrodynamics of this simple modeled system sets clearly in evidence the 
nonlinearity referred to in I, namely, one coming out of the expression for the second flux, cf. 
Eq.(27), and other the higher nonlinear expression for the collision integrals. We also stress 
that in Eq.(26) Bn'(r, t) is also highly nonlinear in the basic thermodynamic variables. 

In that way it has been explicitly shown that nonlinearities are always present in the 
generalized hydrodynamic equations, and we further illustrate the matter first, on the basis 
of the results above, deriving nonlinear equations of evolution for the density and the velocity 
in the fluid under flow, and after that a study of electric current in a system of mobile 
electrons in a doped polar semiconductor, dealt with at the quantum mechanical level. 



III. NONLINEAR EQUATIONS OF EVOLUTION 

Let us consider the equations of evolution for the basic variables of Eqs.(8) to (10), given 
by Eqs.(22), (23) and (24). Deriving in time Eq.(23) and using Eq.(24) it follows that 

§^n{r, t) = - V ■ |l.(r, t) = V ■ V • (r, t) + V ■ - V ■ F,.,(r, t) , (29) 

where for simplicity we have taken the tensorial Maxwell characteristic time of Eq.(26) as a 
scalar (see Appendix A). 
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Equation (29) can be rewritten as 

e„(r, t) |^n(r, t) + ^n(r, t) + V- [n(r, t) v„(r, t)] = 

e„(r, V • V • /Pl(r, t) - e„(r, I„(r, t) ■ VQ-^r, t) , (30) 

after multiplying by © and the definition 

e„(r, t) V • Fe^t(r, t) = V- [n{r, t) v„(r, t)] , (31) 

introducing Va which we call velocity responsible for creating a driven flow in the system. 

Second, in order to close Eq.(30) we need to express the second-order flux in terms of the 
basic variables, n and I„, which after some calculus takes the form 

/PI (r, t) = n(r, t) ks T*{r, t) l^l + n(r, t) [v(r, i)v(r, t)] . (32) 

Using this Eq.(32) in Eq.(30) and, for further simplifying the presentation, let us take 
6„(r, i) as smoothly dependent on r and t, i.e. its gradient and time derivative can be 
neglected, we finally obtain that 

e„(r, t) ^n(r, t) + ^n{r, t) - P„(r, t)VMr, t) = 

-V • [n(r, t) v„(r, t)] + e„(r, i) V • V • (n(r, t) [v(r, Ov(r, t)] ) + ^^(r, t) , (33) 

where G„ contains gradients and time derivatives of ©n(r, t) which, for simplicity we ignore 
in what follows, 

P„(r,i) = ^^;2(r,i)e„(r,i), (34) 

we have used that I„(r, i) = n(r, i) v(r, i), and the second contribution on the right-hand 
side - arising out of the contribution of the so-called convective pressure - can be written as 

V-V - (n(r,i) [v(r,i)v(r,i)]) = (v(r, i) • v)' n(r, i) + 2(v(r, i) • Vn(r, i)) V • v(r, i) 

+ n(r, t) (V ■ v(r, t)) V 2 n(r, t) (w{r, t) v(r, t)) 

+ n(r,t)APl(r,t)®APl(r,t) 



+ 



v(r, t) : Vn(r, t)] ® A'^l (r, i) , (35) 



where A'^l is the velocity-gradient tensor, A'^l = [Vv(r, i)]. 
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Hence, the equation of evolution for the density, Eq.(3), is dependent on the velocity field 
and then it needs be coupled to the equation of evolution for the velocity, which is 

n(r, t) [I + v(r, t) • v] v(r, t) ^ -^V ■ (r, t) - n{v, t) ""^"''I'^^f^ , (36) 

where the pressure tensor field is given by p[^l(r, — n{r,t) kBT*{v,t) 1^1 

Equation (33) is a nonlinear Maxwell-Cattaneo-hke equation with sources. Moreover, 
in the strong condition of neglecting the right-hand side and the second derivative in time 
implying in a movement in conditions such that a;0„ -C 1, we recover the standard Pick's 
diffusion equation 



^n{Y,t)-Vn{r,t) V2n(r,t) =0. 



(37) 



Next let us introduce the effects of advection in the diffusive regime, i.e. the density 
follows the diffusion- advection equation of evolution, as given by Eq.(33), in the condition 

ojQn ^ 1 (thus neglecting the second derivative in time), and for simplicity we illustrate 
the matter in the case of motion restricted to proceed in x-direction; hence we do have that 



d 



d 



— n(a;, t) - V^i^x, t) Q^n{x, t) = -— ^n{x, t) Va{x, t)J + e„(a;, t) 



n{x,t)v^{x,t) (38) 



Moreover, the diffusion coefficient has the standard expression of kinetic theory (in one 

dimension), Vn = v^f^Qn, where Vth is the local and instantaneous thermal velocity 
'mv^f^{x,t) = kBT*{x,t). 

Equation (38) can be reorganized in the following way: 

d d [ 1 d'^v'^ix t) 

— n(a;,t) - Vn{x,t) —n{x,t) = -—\n{x,t)va{x,t)^ +2Qn{x,t)n{x,t)v{x,t) — 



+ 2Qri{x,t)n{x,t) 



dv{x, t) 
dx 



n 2 



+ 2 e„(x, t) v{x, t)^n{x, t)^v{x, t) , (39) 



where 



p„(x,i) = e„(x,i) [v1^{x,t) + v\x,t)\ , 



(40) 



that is, a diffusion coefficient whose standard kinetic theory expression is modified adding 
to the contributions of the thermal velocity the one of the drift velocity. This contribution 
together with the last three on the right-hand side of Eq.(39) have their origin in the second 
divergence of the convection pressure, and it can be noticed that all are nonlinear (quadratic) 
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in v{x,t). The first term on the right of Eq.(39), we recall, is related to the divergence of 
the driving force producing the advective motion [cf. Eq.(31)]. 

Moreover, as already noticed, Eq.(39) needs be solved in conjunction with the equation 
of evolution for the velocity, in this case of motion under flow given by 

d d 
n(x, t)—v(x, t) + n(x, t) v(x, t)—v(x, t) = 
ot ox 

—T*(x,t)n(x,t)—n(x,t)Q^^(x,t) [v(x,t) — Va(x,t)] , (41) 

m ox 

where we used that p{x, t) = n{x, t) ks T*{x, t). 



IV. NONLINEAR QUANTUM TRANSPORT IN SEMICONDUCTORS 

We consider the case of polar semiconductors described by a two-inverted-parabolic bands 
model (in the effective mass approximation and conduction band secondary valleys are 
ignored), where a concentration He of mobile electrons in the conduction band has been 
created by doping. A constant eletric field of intensity Eq in, say, x-direction is applied, 
which accelerates the electrons ("hot" electrons) while there follows a transferring of their 
energy and momentum (in excess of equilibrium) to the phonon field. The sample is in 
contact with a thermal reservoir at temperature Tq, with the phonons being warmed up 
in scattering events involvin g Fr ohlich, deformation potential, and piezoeletric interactions 
with the "hot" electrons P. Il^ (see also Ch. 6 in the book of Ref.[ll]). Scattering by 
impurities is neglected in comparison with the one due to lattice vibrations. Moreover, in 
these polar semiconductors of the different types of electron-phonon interaction we keep only 
the predominat Frohlich (polar) interaction. 

The electric field E creates a (uniform in space) current, that is, a flux of charged particles, 
and then we have a kind of hydrodynamic motion of the like of the one analysed in the 
previous sections, except that, as noticed, the hydrodynamic variables are uniform in space 
(independent of position r), and the thermal bath is played by the phonons. We resort to 
a first-order hydrodynamics, as in section II, namely, we consider the carriers' energy, the 
density of particles and the flux of particles. The latter multiplied by the effective mass of 
the conduction band electrons is the linear momentum, and while the density of particles 
remains constant in time, as noticed the energy and the momentum change in time. The 
Hamiltonian is composed of the energies of the electrons in the conduction band (treated as 
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usual in the random-phase approximation thus including Coulomb interaction in a mean- 
field approximation) and that of the polar (longitudinal optical) phonons, we call Hq this 
part, plus the Prohlich interaction of electrons and polar phonons and the interaction of the 
electrons with the electric field, we call H' this part. Moreover, differently in other aspect 
with the case of section II we need here to use a quantum approach. 

Taking as basic dynamical variables Hg, Hlq, P (energy of electrons He and polar 
phonons Hlo, number of electrons, and linear momentum of the electrons) the auxiliary 
("instantaneously frozen") statistical operator is then given by 

p{t, 0) = exp [-(f){t) - m [He - i^{t) N - v{t) ■ P] - f3Lo{t) Hlo} Pr , (42) 

where, we recall, 0(i) ensures the normalization condition, and we have introduced the 
nonequilibrium thermodynamic variables 

[P{t) = l/kBT*{t),-P{t)fx{t),-P{t)v{t),PLo{t) = 1/kBTloit)], (43) 

interpreted as nonequilibrium quasitemperature, T*, quasi-chemical potential, fi, and drift 
velocity, v, of the electrons, and quasitemperature of the polar phonons, TIq. Finally, pr is 
the statistical operator of the reservoir at temperature Tq. 
We call the basic macrovariables 

[Ee{t),N,P{t),ELo{t)], (44) 

where is constant in time, and we recall, the density is indicated by n, and the nonequi- 
librium thermodynamic equations of state are, after direct calculation, given by 

^Ee{t) = ^nknT^t) + ^nm*v^{t) , (45) 
ip(i)=nm*v(i), (46) 

1 1 r 1-1 

- ELo{t) = cuo ex.p{PLo{t)hu;o}-l , (47) 

V Vcell L 

where, 1. In the conditions of concentration and values of quasitemperature involved in 
the applications, the nonequilibrium time-dependent Fermi-Dirac-like distribution of the 
internally thermalized electrons in band states in the effective mass approximation can be 
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approximated by a nonequilibrium time- dependent Boltzmann-Maxwell one , and then 
there follows the first term on the right of Eq.(45), namely, equipartition at each time t in 
terms of the quasitemperature T*(t) of the "hot" electrons of effective mass m*, and the 
other term is evidently the kinetic energy of drift (current). 2. The polar phonons have 
been treated in a dispersionless approximation (Einstein model), with unique frequency Uq, 
Vceii is the volume of the unit cell in the crystal, and V is the volume of the sample. 
On the other hand the equations of evolution (see Appendix B) are 



for the energy of the electrons, 



d_ 

di 



P{t) = enE- |Jp(t)| , 



for the electron linear momentum in the direction of the applied field, 

J^ELoit) = \Jpit)\- \Janit)\ , 



(48) 



(49) 



(50) 



for the energy of the polar phonons; volume V has been taken as 1. 

In Eq.(48) the first term on the right accounts for the pumping of energy on the carrier 
system because of the presence of the electric field, while the second represents the rate of 
excess energy transferred to the polar phonons. In Eq.(49) the first term on the right is the 
force produced by the presence of the electric field, and the second the rate of momentum 
transferred to the lattice. In Eq.(50) the first term on the right is the gain of energy pumped 
on the phonons by the nonequilibrated ("hot") carriers, with the second being the transfer 
- via anharmonic processes - of such energy to the acoustic phonons acting as a thermal 
bath (finally, from the latter there follows a transfer to the thermal reservoir). As noticed, 
their expressions a re g iven in the Appendix B with the calculations being performed in the 



Markovian 



m 



25 



26|. We write down here the two collision integrals of Eqs.(48) and (49) 
to illustrate the highly nonlinear dependence of them on the nonequilibrium thermodynamic 
variables, (3{t) and v(t), which, on the other hand, are related to the basic variables through 
the also nonlinear equations of state, viz. Eqs.(45) and (46): 



Mt) I = y^^^ e-" n uoA^- (z/q + 1) A2 



(51) 



Jp(t) I = Ap 1/^/2 e"^ n z/o A3 - (z/o + 1) ^4 



(52) 
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where 

y{t) = m ^^0 ; x{t) = v\t)/kB T*{t) , (53) 



A, 



2hujQ 



vrm' 



1/2 



A, 



20F 



(54) 



^0 



exp {/?Lo(i) ^t^o} - 1 



-1 



(55) 



Ai{t) 



K,{y/2) + Y, 



(4xy)" 
{ {2n + l)\ 



■IT 



Ko{y/2) 



(56) 



n=l 



K,{y/2) 



(57) 



1 , r 

A,{t) = {Axyy/'— Ko{y/2)-K^{y/2) 
3 y I 



'n— 1 



n=2 



{2n + l)\ 



dy 



n-l 



iro(l//2)-i^i(z//2) 



(58) 



1 , ev/^ r 

A,it) = e-y{AxyY/^— K,{y /2) + K,{y /2) 
3 y L 



n=2 



2n 



(2n + 1)! 



(4x?/)' 



r 



'71—1 



n-l 



;y/2 



Koiy/2) + K,{y/2) 



■ (59) 



In these expressions - and note that we have omitted on the right-hand side to exphcitly 
indicate the dependence on time of the different quantities involved - stands for Bessel 
function of order n [l^ and Sq is the Frohhch field intensity in the polar interaction isj ]. 

It is worth noticing that in quantities ^3 and A4 [cf. Eqs.(58) and (59)] it can be put in 
evidence the quantity x^/^ which is proportional to v, and therefore we can formally write 
Eq.(49) as 

d 



/ N e ^ v(t) 
-v{t) = —Eo - 
dt m* Tp [t) 



(60) 



after using Eq.(46), being an equation of the Newton-Langevin type but where Tp(t) plays the 
role of a relaxation time of the velocity (momentum) which is depending on time through 
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its dependence on the nonequilibrium thermodynamic state of the system, i.e. depends 
on a highly nonhnear way on PLoit), and the velocity v{t). Equation (60) can be 

alternatively written in the form of the integral equation 

^W = ^^or(t), (61) 

fit 



where 



POO 

r,(t)=exp{-^(t)} / dt' expi^Pit')} , (62) 

^0 



with 

ft 



m= / dt'T-\t'), (63) 



after taking into account the initial condition f (0) = 0. Moreover, once the current density 
is given by 



I{t) = env{t), (64) 



using Eq.(61) we can write 



I{t) = a{t)Eo, (65) 



introducing a time-dependent (on the evolution of the nonequilibrium thermodynamic state 
of the system) Drude-type conductivity 

a{t) = ^rM- (66) 

m* 

Next we illustrate numerically the matter considering the polar semiconductor GaAs and 
the strong polar ones and large gap GaN (of present interest for its use in blue diodes and 



lasers 1J|). Figure 1 shows the increase, in the steady state (which follows after a transient 
of a few picoseconds), of the electrons' quasitemperature with the electric field strength in 
the case of doped {n ~ lO^^cm^^) GaAs at a reservoir temperature To = 3007^. The dots are 
the result of a Nonequilibrium Molecular Dynamics (Monte Carlo-style) simulation (from 
Ref. [15]), where we can see a good agreement between both types of approaches. 

In Figure 2 are diplayed the results of the calculations of the electrons' mobility (left 
ordinate) in the steady state of doped (n ~ 5 x lO^^cm"^) GaAs, as well as the momentum 
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relaxation time (right ordinate), in the Ohmic regime, as a function of the reservoir tem- 
perature: It can be noticed a good agreement with the experimental data taken from three 
sources, namely, Refs.[16], [17] and [18]. 

Going over the particular case of the large-gap strongly polar GaN (in the cubic, i.e. 
zincblende, phase with n ~ IQ^^cm^^ , and Tq = SOOfC), we can see in Figure 3 the increase 
of the electron- drift velocity with the electric field intensity, and a good agreement can 
be noticed in a comparison with the NMD-Monte Carlo simulation, the dots, taken from 
Ref.[19]. 

Finally, in Figure 4 we can accompany the evolution in time of the electron-drift velocity 
in doped {n ^ IQ^'^cm"^) cubic GaN, where in the horizontal axis we do have the travelled 
distance at any time t, i.e. ft, and in the presence of an electric field intensity of ?)QKV / cm. 
We can notice, first, that the transient time is roughly of the order of 300 femtoseconds 
and it can be observed the presence of a so-called "velocity overshoot" at roughly 100 
femtoseconds. The agreement with the NDM- Monte Carlo simulation, of Ref.[20], is very 
good. The calculations and figures are from Ref.[21]. 



V. CONCLUDING REMARKS 



As noticed in the Introduction, in the previous article in this issue it was described the 
construction of a generalized higher-order-nonlinear hydrodynamics based on a nonequi- 
librium ensemble formalism. This means a mechano-statistical foundation in terms of a 
generalized nonequilibrium grand-canonical ensemble, in that way establishing a kind of 
unification of the microscopic dynamics and a mesoscopic hydrodynamics. In this follow up 
article we present some illustrations on the working of the formalism. In Section II we have 
considered the case of a system composed of two ideal fluids in interaction between them. It 
is used a truncated hydrodynamic description in which are included solely the density of en- 
ergy, the density of particles, and the first flux (current) of particles. As discussed elsewhere 
this implies on restrictions on the characteristics of the hydrodynamic motion, basically the 
case of those restricted to long to intermediate wavelenghs jslls^, i.e. a hydrodynamics of 
first order. The continuity equations are derived [cf. Eqs.(22) to (24)] and to close them, in 
this truncated description, one needs to express the second-order fiux in Eq.(24) in terms 
of the basic variables. In doing this we find a first source of nonlinearity as can be seen in 
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Eq.(27). Moreover, the collision integrals on the right of these equations of evolution are 
highly nonlinear as shown, for the case of the one in the equation of evolution for the first 
flux, in Appendix A in the model of two ideal fluids in interaction, and in Appendix B for 
the case of a fluid of conduction electrons in doped semiconductors. In the first case the 
collision integral is composed of several contributions: One dependent on the^ gradient of 
the density, as it should appear in a mesoscopic hydrodynamics of this type 23]. Other ac- 
counting for nonlocal effects, and a local in space third one which is of the type of Maxwell's 
contribution 0,0,0]; cf. Eq.(A6). It is of a tensor ial (rank 2) character and highly nonlinear 
in the nonequilibrium thermodynamic variables, and then dependent on the position and 
time. 

Section III deals with the motion under flow (advective motion) for the case of the 
description of Section II. The equations of evolution for the density of particles and for the 
field of velocity - which are coupled - are obtained. In that way are derived, respectively, a 
generalized Maxwell-Cattaneo-type equation and a generalized Navier-Stokes-like equation. 
The nonlinearities are present, and the generalized diffusion coefficient is dependent on the 
local and instantaneous nonequilibrium thermodynamic state of the system. 

Finally, in Section IV it is considered quantum transport in doped polar semiconductors, 
that is we deal with the homogenous current of the mobile electrons. In this case it is obtained 
in the steady state - that sets in after a very short (nanosecond scale) transient - when 
under the action of a constant electric field (that is, the space dependence carried on in the 
previous Sections is not present). It is derived the dependence of the nonequilibrium electron 
temperature (carrier's quasitemperature) on the electric field intensity as well as of the drift 
velocity (the current is proportional to it) in several cases: Comparison with computational 
modeling Monte Carlo calculations and with experimental data is done following a very good 
agreement. 

Acknowledgment s 

We acknowledge financial support from Sao Paulo State Research Foundation (FAPESP). 
ARV and RL are Brazil National Research Council (CNPq) research fellows. 



15 



APPENDIX A: THE COLLISION INTEGRAL 



The collision integral present on the right-hand side of Eq.(17), resulting from the colli- 
sions of the two types of particles via the potential V{\rj — R^|) [cf. Eq.(l)] is given in the 



Markovian approximation 

Jn 



25 



26[ by the expression 



r,t)~ dt'e<''-''> j dT [H'{T\t' -t)o,{H'{T),UT\r)}} p{t,0) X pn, (Al) 

where H'{t' — t)o indicates evolution under the dynamics generator by Hq - i.e. in the 
interaction representation -, and we use for the statistical distribution of the reservoir, p/j, 
a canonical one with temperature (3^. It can be noticed that the approximate (Markovian) 
scattering integral of Eq.(Al) is quadratic in the interaction strength and corresponds in 
the classical limit to the Golden Rule of Quantum Mechanics, and we recall that {...,...} 
stands for Poisson bracket. 

The lengthy but straightforward calculation of Eq.(l) provides several contributions to 
this scattering integral: one involving the gradient of the concentration (leading to a correc- 
tion (renormalization) of the diffusion coefficient due to the collisions), other involving the 
gradient of the quasitemperature (implying in a cross-effect between thermal and material 
motion). Another one is related to introduce correlation effects in space of the coUisional 
processes (i.e. nonlocal effects). Finally, the most relevant contribution is given by 

Nr /tt .^„^\V{q)\^ 



Jn(r,t) 



mrrir \l 2 



(A2) 



where 

A(r,q) = n(r,t) 



27rm 



3/2 



d^p p exp 



[p -mv(r,t)]^ - 



MI3r 



q- p 



(A3) 



In Eqs.(A2) and (A3), V(q) is the Fourier transform of the interaction potential between the 
two types of particles, is the reduced mass, m"""^ = m^^ + M~^, and is the number 
of particles in the thermal bath. 

Introducing a shift in coordinate p, namely 



P + b 



with 



b(r, t) = /5 m v(r, t) 



(A4) 



(A5) 
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and performing the calculations we arrive at the result that 



(A6) 



where G is a rank- two tensor with dimensions of time, playing the role of a tensor ial Maxwell- 
relaxation time , with its inverse given by 



mr,t)] 



-1 nR/3Rm + M /ttM/Sr 



m m 



1/2 



m/5(r, t) 



3/2 



where 



m/3(r, t) + M(3r 

i^[2](r,t) = l^|V(g)|2 ^ mM/3ij/?(r,t) (q-v(r,t))2 



K^'\v,t), (A7) 



exp 



2[m(3{r,t) + Mf3R] 



(A8) 



with ur being the density of particles in the thermal bath. 

We stress that this tensorial Maxwell-relaxation time depends on position and time 
through its dependence on the basic nonequilibrium thermodynamic variables /3(r,t) and 
v(r, t): recalling that In(r, t) = n{r, t) v(r, t), this time-relaxation-type contribution is not a 
linear one, but a highly nonlinear in v(r, t). In fact, performing the integration in Eq.(A8) 
we find that Kij = for i ^ j, and there survive the diagonal terms 



K,,{r, t) = Kyyiv, t) = m $ (^^(r, t) v\v, t)) - ]^K,,{j, t) 



^ ~' ^ 2/ n1 "'^^^ / 5/ \ 2/ 



where 



/3(r,t) 



M(3R(3{v,t) 
m(3ir,t) + M(3R' 



k=l 



m 



{2k-l){k-l)\ \ 2 



-/?(r,t) v\v,t) 



90 



|V(q) 



(A9) 



(AlO) 



(All) 



(A12) 



(A13) 



with go being a cut-off limit of the order of the inverse of the mean distance between particles 
of the system and the bath, 7 is the incomplete Gamma function I4I, and we can also 
alternatively write 



i^..(r,t)=e(0)$^(-l)'= 



(fc + l) k\\2 



(A14) 
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It may be noticed that we can write 



l[2] 



e„(r,t) 



+ 



o[2] 



©n (r,t) 



(A15) 



where 



3^ = iT.{Kl(M)]-}. 



(A16) 



and the last term is then traceless, and accounting for anisotropic effects. Furthermore, 
using Eq.(AlO) there foUows that 

= - [KUr. t) + K,y{v, t) + i^,,(r, t)) = -m $ (-/3(r, t) ^^(r, t) j . (A17) 

Conserving only this contribution in Eq.(A14), hnearizing in v(r, all expressions, say, the 
case of a low density flux, assuming a good thermal contact between system and reservoir, 
such that /9(r, t) ~ (3r, and small amplitude movement, namely n(r, t) ~ n, we are left with 
a typical Maxwell-type contribution 



J„(r,t) ~ -ln{r,t)/Qn, 



(A18) 



where 



1 



\l m{m + M) 



(A19) 



APPENDIX B: ENERGY AND MOMENTUM RELAXATION, EQS.(47) AND 
(48) 



The collision integrals are in this case given by 



J?\i)-i^y l^dt'e<''-'^Tr[[H'{t'U p(i,0)}, (Bl) 



where Ai = for energy and ^2 = P for momentum, and H' is Prohlich electron-phonon 
interaction; subindex nought indicates time dependence in interaction representation. These 
operators are 



(B2) 
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kq 

P = 5^^kc[c„ (B4) 

k 

where (c|^) are annihilation (creation) operators of electrons in band state k, and the 
effective mass approximation is taken; (^^) ^"^^ annihilation (creation) operator of LO 
phonons in mode q, and the matrix element in Frohlich interaction is 

Cq = -i^, (B5) 

a being Frohlich coupling constant. 

After some calculation one arrives at the results that 

kq L 

- m (1 - - - ^^,) |, (B6) 

kq I 

- ^.{t) m (1 - A+cW) '^(^.+. - - ^u;,) |, (B7) 
where, in a nondegenerate-like limit acceptable in the usual experimental conditions 

at) = A{t) exp |-/5(t)^ (k - kj^l , (B8) 

a shifted Maxwell-Boltzmann-like instantaneous distribution, with ^k^ = mv, P = nv, 
A{t) ensures the normalization to N , ujq is the LO phonons dispersionless frequency (Einstein 
model), and [cf. Eq.(54)] 

{t) = [ exp {(5lo {t) huJo]-l\ \ (B9) 

It can be noticed that Eqs.(B6) and (B7) are of the form of the Golden Rule of Quantum Me- 
chanics averaged over the noequilibrium ensemble. Performing the integrations one arrives 
to the results presented in Eqs.(47) and (48). 
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FIGURE CAPTIONS 




FIG. 1: The electron quasitemperature vs. the electric field intensity, in the steady state of n- 
GaAs, comparing the results of the NESEF-based kinetic theory with a Monte Carlo simulation 
(full circles from ref.[15]). 
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FIG. 2: Electron mobility (and momentum relaxation time) in n-GaAs for different values of the 

reservoir temperature, comparing tlie results of the NESEF-based kinetic theory with experimental 
data (up triangle from Ref.[14], x from Ref.[17], empty circles from Ref.[18]). 
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FIG. 3: Electron-drift velocity vs. electric field intensity in n-GaN, comparing the results of the 

NESEF-bascd kinetic theory with a Monte Carlo simulation (Full circle from Ref.[19]), with Tq = 
300i^. 
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FIG. 4: Evolution of the electron-drift velocity in 77,-GaN in terms of the travelled distance, com- 
paring the results of the NESEF-bascd kinetic theory with a Monte Carlo simulation (Full circles 
from Ref.[20]): zincblende n-GaN for an electric field intensity 30kV/cm and Tq = 300K. 
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